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ABSTRACT 






A computer simulation of the Brayton cycle is used to develop a 
design procedure with respect to minimizing volume and weight for the 
counterflow plate-fin recuperator. Based on the Orenda OT-4 600 horse- 
power gas turbine, recuperator performance and dimensional character- 
istics are presented for an idealized equilateral plate-fin recuperator 
core matrix. The effects of plate spacings varying between 0.1 inches 
and 0.5 inches on recuperator performance characteristics are presented 
over a wide range of core frontal areas. Specific trends toward minimum 
volume and minimum weight plate-fin recuperator cores are discussed in 



detail . 
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NOMENCLATURE 



English Letter Symbols 



A 

A. 



fr 



a 

b 



C 

F 

G 

g 

h 

HP 

i 

k 

L 

m 

P 

P 

w 



heat transfer area 
polynomial coefficient 
free flow area 
frontal area 
plate thickness 
plate spacing 
recuperator width 
specific heat 
heat capacity rate 
force 

mass flow rate per unit area 

acceleration of gravity 

convection heat transfer 
coefficient 

horsepower 

specific enthalpy 

thermal conductivity 

flow length 

mass flow rate 

pressure 

wetted perimeter 

heat transfer rate 

lower heating value of fuel 



sq ft 



sq ft 
sq ft 
ft 
ft 
ft 

BTU/ (lbm deg F) 

BTU/ (hr deg F) 
lbf 

lbm/(hr sq ft) 

(lbm f t) / ( lbf sq sec) 

BTU/ (hr sq ft deg F) 
hp 

BTU/ lbm 

BTU/(hr sq ft deg F)/ft 
ft 

lbm/ ft 
lbf /sq ft 
ft 

BTU/hr 
BTU /lbm 
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English Letter Symbols (continued) 



r h 


hydraulic radius 


ft 


s 


fin spacing 


ft 


S 


specific entropy 


BTU/(lbm deg R) 


SFC 


specific fuel consumption 


lbm/(hp hr) 


T 


tempe rature 


deg F, deg R 


U 


overall heat transfer 
coefficient 


BTU/ (hr sq ft deg F) 


U 


velocity 


ft/sec 


V 


volume 


cu ft 


Vs 


solid volume 


cu ft 


V 


specific volume 


cu ft/lbm 


w b 


available shaft work 


BTU/lbm 


W 

c 


compressor work 


BTU/lbm 


Wt 


mass 


lbm 


w 


fin thickness parameter 


ft 


Greek Letter 


Symbols 




oC 


ratio of total transfer area 
of one side of exchanger to 
total exchanger volume 


sq ft/cu ft 




compactness 


sq ft/cu ft 


6 


fin thickness 


ft 


c 


regenerator effectiveness 




& 


fin angle 


deg 


Go 


temperature difference 


deg F 




fluid viscosity 


lbf hr/sq ft 


e 


density 


lbm/cu ft 




ratio of free flow area to 
frontal area 


sq ft/sq ft 
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Efficiency Terms 





combustor efficiency 




compressor efficiency 




fin effectiveness 


’V 


thermal efficiency 




turbine efficiency 




surface effectiveness 



Dimensionless Parameters 



$ 


Fanning friction factor 


J 


Colburn-j factor 


K 

c 


contraction loss coefficient 


K 

e 


expansion loss coefficient 


m 


fin heat transfer parameter 


NTU 


Number of Transfer Units 


Prel 


relative pressure 


P 


relative pressure drop 


Pr 


Prandtl number 


Re 


Reynolds number 


r f/a 


fuel to air ratio 


St 


Stanton number 
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I. INTRODUCTION 



The heat exchanger, a device which effects the transfer of heat from 
one fluid to another, is found whenever it is necessary to transfer 
energy without mixing the fluids involved. It is felt by many that our 
future as a power rich nation may depend partially on the heat exchanger 
engineer, as the most promising sources of power often present the most 
difficult applications. One such power source, the gas turbine, has a 
number of inherent advantages over the common gas and diesel engines; 
among them are compactness, low weight, ease of maintenance, and multi - 
fuel capabilities . 

With the addition of a heat exchanger of reasonable size, the ther- 
mal efficiency of a gas turbine unit may be increased to about the same 
as, or slightly higher, than the gas engine [7].* The added weight of 
the exchanger, however, could well negate one of the primary advantages 
of the gas turbine. Furthermore, the relative rise in production cost 
due to the addition of a heat exchanger ought to be considered in any 
comparison, so as to keep the gas turbine on a specific cost level with 
the diesel engine. 

Heat exchangers, however, are inherently bulky, heavy, and expen- 
sive. To overcome this problem, the concept of the compact heat ex- 
changer has come into being. Not necessarily small, such an exchanger 
incorporates a heat transfer surface of high area density (leading to a 
high effectiveness for given weight and volume limitations) and, generally, 

* 

Numbers in brackets designate references listed in Bibliography. 
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may be manufactured at a much lower cost per square foot of heat trans- 
fer area than the not so compact heat exchanger [ 9 ] . 

Although standard lines of compact heat exchangers have been adopted 
by many manufacturers, such exchangers may not be purchased as off-the- 
shelf items, as a proper design is needed for nearly every requirement* 
Such a design procedure usually considers the factors of heat transfer 
requirements, pressure drop characteristics, size, and cost, with an 
emphasis on the balance between the gain due to heat recovery and the 
loss due to pressure drop contributed by the addition of the exchanger* 

There exist, however, numerous applications for which this much 
used design philosophy is not sufficient* Often the designer is faced 
with stringent volume or weight limitations. The technique of deter- 
mining an optimized counterflow heat exchanger with respect to volumetric 
and weight considerations is developed in this study. 

Fundamental relationships, assumptions, and design points upon 
which the investigation is based are presented in Sections II through 
IV; the computer program used to simulate the regenerated gas turbine 
cycle is discussed in Appendix B. 
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II. THE IDEAL REGENERATED BRAYTON CYCLE 



The gas turbine plant, to a large extent, must have all its units 

designed with reference to the cycle on which it is to operate. The 

designer, therefore, must have a knowledge of the particular cycle to 
be considered. The most prevalent type of gas turbine in use today, 
the simple open cycle, has achieved great success due to its versatility 
[6]. With the addition of regeneration the efficiency of this cycle may 
be greatly increased. 

The Bray ton cycle, the basic thermodynamic cycle of the gas turbine, 
(refer to Figure 1) is characterized by taking air, as the working fluid, 
from the atmosphere into a compressor where part of the turbine power is 

used to raise the pressure before the air enters the cold side of a re- 

cuperator. The recuperator reduces the heat input to the cycle by trans- 
ferring heat from the hot turbine exhaust gases to the cooler air. A 
further temperature increase occurs as fuel is burned directly in the 
air flow (forming gaseous products) through the combustor. The gas then 
enters a turbine, where it expands to a lower pressure while releasing 
energy available for work. Upon leaving the turbine stage, the gas 
enters the recuperator hot side and transfers much of its otherwise waste 
heat to the incoming air, before being exhausted into the atmosphere. 

The particular cycle used in this heat exchanger investigation 
is based on the Orenda OT-4 gas turbine, a 600 horsepower unit engine- 
ered for the U. S. Navy [10]. The design specifications for the Orenda 
unit are: 
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HP 


net power output 


600 hp 


V 


overall thermal efficiency 


34% 


pi 


compressor inlet pressure 


2120 lbf/sq ft 


Tl 


compressor inlet temperature 


560 deg R 


T4 


turbine inlet temperature 


2200 deg R 


P6 


exchanger exhaust pressure 


2120 lbf/sq ft 


P2/P1 


compressor pressure ratio 


4.0 


P3 - P4 


combustor relative pressure drop 


0.02 


P3 








compressor efficiency 


86% 


Q 


combustor efficiency 


98% 




turbine efficiency 


86% 


T f 


fuel inlet temperature 


530 deg R 




fuel lower heating value 


19000 BTU/lbm 



Other parameters which are usually specified for the cycle, such as the 
mass flow rate, regenerator effectiveness, and exchanger pressure drops, 
are not presented; they are, in fact, considered as variables in the 
optimization procedure to be discussed. 

A more detailed discussion of the Brayton cycle is outlined in 
Appendix A. 
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Temperature, T 
Enthalpy, i 




Entropy, s 







FIGURE 1 THE IDEAL REGENERATED BRAYTON CYCLE 
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III. RECUPERATOR SURFACE AND CORE GEOMETRY 



For all applications in which the turbine inlet temperature is 
limited and a low to moderate pressure ratio is to be used, the most 
efficient method of improving the economy of the gas turbine involves 
the use of a heat exchanger, in which hot turbine exhaust gases give up 
heat to the relatively cool air from the compressor. Two basic heat ex- 
changer models, the recuperator and regenerator, are in common use. The 
term recuperator (which is considered in this study) is usually applied 
to a continuous flow method of heat exchange, while the term regenerator 
specifies a periodic flow device, usually involving a moving matrix. 

A. COMPACT HEAT EXCHANGERS 

Considerations involving heat transfer characteristics of the gas 
turbine recuperative unit have led to the development of the compact 
heat exchanger. The term compact does not refer directly to a small 
bulk or weight, but to the heat transfer surface density of the ex- 
changer. The most useful measure of this density is the compactness, 

(2^ , whe re 



Somewhat arbitrarily, the working limits of compactness for a simple 
counterflow recuperator may be specified as 1000 sq ft/cu ft as the 
maximum working limit, and 200 sq ft/cu ft as the minimum limit before 
leaving the realm of compactness [9]. 

An understanding of the geometric parameters, presented in Figure 
2, and subsequent relations most often used is essential to any compact 
heat exchanger design procedure. 




= total transfer area of one side of exchanger 
volume between the plates of the same side 
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Gases 

FIGURE 2 

COMPACT COUNTERFLOW RECUPERATOR CORE 



a 

b 



r h 

£ 



A 

A 



c 




L 

V 

V 

s 



oC 



plate thickness 

plate spacing 

recuperator width 

flow passage hydraulic radius 

compactness 

heat transfer area of one side 
free flow area of one side 
frontal area 
flow length 
volume 

solid core volume 

ratio of total heat transfer 
area of one side to the total 
exchanger volume 

ratio of free flow area of one 
side to frontal area 



ft 

ft 

ft 

ft 

sq ft/cu ft 
sq ft 
sq ft 
sq ft 
ft 

cu ft 
cu ft 

sq ft/cu ft 
dimensionless 
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The following equations give the relations between surface and 
core factors for one side of the counterflow plate-fin recuperator* 
Subscript 1 refers to any one side, and 2 refers to the other side. 
Factors without a subscript are common to both sides. 



a - _A± - At 
^ ~ V, ~ bxLc 



3.1 



l hi 




L Ac 

A 



1 

- 1 






3.2 



<A 



Ai 

V 



bx fox 

bx * b z *2.0. 



3.3 



<*1 



Ac i - A-l fW l _ / ■p — Pnx bx (bx. 

A; r 1 ' n ' 1 bi*b2.*£.a_ 



A 



CL 




Ai fwi 

L 



3.5 



Al — 



V 






(A-fr) (A-cl 




3.7 



Among the factors common to both sides is A^, the total recuperator 
frontal area. The preceding relations are sometimes written on the 
basis that, for the counterflow recuperator, the hot and cold side 
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frontal areas are defined separately such that their sum equals A^. 

This investigation requires no such separation of terms and is based 
on the interpretation of frontal area as presented in Figure 2. 

B. PLATE - FIN SURFACES 
1 . General Theory 

Although other matrix forms are used, the plate-fin recuperator 
matrix is the most common [6]. Not only does this type of structure lend 
itself to light weight, extremely compact designs, but the plate spacing 
on each side may be easily optimized independently, allowing for design 
flexibility. 

An understanding of the specific geometric parameters of the 



plate-fin matrix (Figure 3) and their relations is necessary to the 
design procedure. 




FIGURE 3 

PLATE - FIN GEOMETRY 



21 



a 


plate thickness 


ft 


b 


plate spacing 


ft 




fin thickness 


ft 


S 


fin spacing 


ft 


A 


heat transfer area 


sq ft 


A 

c 


free flow area 


sq ft 


P 

w 


wetted perimeter 


ft 


L 


recuperator length 


ft 


w 


fin thickness parameter 


ft 



Cj — 2-Cb - w") 

Ton© 



Ac- “2" 



p ~- zl * s - a(a^- e )+ £> - ^(dr e + i) 3.10 



A 




3.11 



From these relationships, the values for two very important parameters 
may be found. 




- (b-w)[2C 



i 




3.12 
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$ = -^[( a b *~J 



bS 




b(c^e + 1) 
b 



3.13 



Since the recuperator matrix studied involves a fin angle, 0 , 



of sixty degrees (resulting in an equilateral fin matrix) a further 



simplification may be made. 



w = 



3.14 



, _ (b-20 



3.15 






j o 

b 



3.16 



Table I represents values of these parameters for the plate spacings 
investigated. Although the spacings greater than 0.4 inches fall out 
of the previously defined compactness region (200 ^ (b — 1000), they 

are included in the study since: (1) this region is arbitrary and (2) 
recuperator trends in weight and volume are more easily observed when 
large and small spacings are compared. 

2 . Specific Parameters 

In order to keep the investigation on a common basis, all cores 
analyzed are assumed to be constructed of stainless steel, with con- 
stant properties, and uniform plate and fin thicknesses. The constant 
parameters used in the analysis are: 



a 



plate thickness 



0.0008333 ft 
(0.01 in) 



fin thickness 



0.0004166 ft 
(0.005 in) 
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k 



f 



Pf 



thermal conductivity of 
stainless steel plates 
and fins 

density of stainless 
steel plates and fins 



12 BTU/ (hr sq 
ft deg F)/ ft 

484 lbm/cu ft 
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TABLE I 



EQUILATERAL PLATE - FIN SPACINGS 



Plate 


Compactness, P> 


Hydraul i c 


Spacing, b 


\ 


Rad i us, r h 


(In) 


(sq ft/cu ft) 


(ft) 


0.1 


720.00 


0.0012500 


0.125 


576.00 


0.0015972 


0.15 


480.00 


0.0019444 


0.175 


411.43 


0.0022917 


0.2 


360.00 


0.0026389 


0.25 


288.00 


0.0033333 


0.3 


240.00 


0.0040278 


0.35 


205.71 


0.0047222 


0.4 


180.00 


0.0054167 


0.5 


144.00 


0.0068056 




* 

The hydraulic radius. 


r. , is based 
n 




on a fin th i ckness, 6 , 


of 0.005 in. 



25 



IV. HEAT TRANSFER AND FLOW -FRICTION DESIGN THEORY 



The design of a recuperator involves a consideration of both the 
heat transfer rates between the fluids and the pumping power expended to 
overcome fluid friction in the exchanger [8]. Obviously, one of the most 
desirable characteristics of a compact heat exchanger is that it yield a 
maximum heat exchange for a minimum pressure drop. 

A. HEAT TRANSFER CHARACTERISTICS 

Elementary heat transfer principles show that the cool air and hot 
gas flows should be in opposite directions ; otherwise, the two streams 
would approach their mean temperature only. Subsequently, many compact 
heat exchangers involve a counterflow design, with plate-fin construc- 
tion for added efficiency [5]. 

An overall heat transfer rate equation which includes convective 
and conductive heat mechanisms responsible for the heat transfer from 
the hot gases to the cool air may be written for a differential element 
as 



where U is an overall heat transfer coefficient, dA represents an 
element of heat transfer area, and dq is a differential heat flux. 

A physical understanding of the heat transfer performance of the heat 
exchanger may be attained with the use of equation 4.1 and a number of 
parameters to be defined. 

1. Fin Effectiveness 




4.1 



The straight recuperator fin of constant thickness & and 
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length Jl is represented in Figure 4. Such a fin extends from the 

base plate, at temperature T^, through the moving fluid, at temperature 

T , to a second plate. As there is no heat flux through the fin 

oo 

midway between the plates (due to the symmetry involved), a heat trans- 
fer analysis need only consider a fin of effective length Si /2, insul- 
ated at the tip (Figure 4) . An energy balance over an element of this 
fin [2] results in the equation: 

q = k m A f © 0 Tanh(m 4.2 



k = 
m Ct 



thermal conductivity of fin material 




= a dimensionless parameter for thin fins 



^ & C 1) 



cross section area of fin per unit width 







- T 



<y> 




FIGURE 4 

STRAIGHT FIN OF UNIFORM THICKNESS 
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In the order to express the heat exchanging capacity of an extended 
surface it is useful to define the term fin effectiveness as 



it 

ns. 



fin effectiveness 

actual heat transferred from fin 

heat transferred if entire fin were at temperature T 



0 



k m A c Oo TanhQm JL /2 ) 
2 L h ( A /2) 



= 

= (Tanh(m JL /2)) (m JL /2) 



-1 



4.3 

4.4 



2 . Surface Effectiveness 

Since the plate-fin geometry involves an array of fins, 
it is useful to define a total surface temperature ef fectiveness , , 
which gives a measure of performance of the entire array. Defining 
M A M as the total heat transfer area of the combined fin and plate 
surfaces , 



_ actual heat transfer from A 
heat transfer if A were at T 



0 



^ (A - A f )h(T 0 - TJ+ A f h(T Q - Too ) 



4.5 



4.6 



A h (T q - T <p ) 



*lc 






3 . Overall Heat Transfer Coefficient 

The overall conductance U, as defined previously, may 
now be calculated as 



Ut - [ 



(Ah/AORo.h'Hv 



CAw/A 0 k 



1 



4.7 



where U is based on a unit of cold side total surface [2]. In the 
c 
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usual gas-to-gas recuperator, the second term in equation 4.7 (the wall 
resistance component) may be neglected relative to the fluid resistances. 




overall surface effectiveness for hot and 
cold sides 



A 



w 



average wall area 



A, 



h 



total heat transfer area of hot and cold sides 



A 



c 



h, 



h 



convection heat transfer coefficients of 
hot and cold sides 



h 



c 



a 



plate thickness 



k 



thermal conductivity of plate material 



4. Counterflow Exchanger Heat Transfer Effectiveness 

As previously stated, the most common and most effective flow 
arrangement for a two-fluid heat exchanger is the counterflow design. 
Referring to Figure 5, an energy balance leads to a second equation for 
the amount of heat transferred between the two fluids. 



Integration over the length of the recuperator leads to the equation 



,ia 'T'h.ooO — CLp C^c.oo t 4 . 9 



It is useful to define the term "capacity rate" as the product of the 




4.8 



specific heat of one fluid and the mass flow rate of that fluid, 

referring to the hot gas, and C referring to the cooler air. 

c 



- Cv\ - Ce pot - ^ 



4.10 
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T h, in 



T 

c,out 




h ,out 

T . 
c, in 



FIGURE 5 

TEMPERATURE VARIATIONS IN COUNTERFLOW HEAT EXCHANGER 

The heat transfer performance of the recuperator may now be deter- 
mined by the heat exchanger effectiveness, £ , defined as the ratio of 
the actual heat transfer to the thermodynamic maximum heat transfer 
rate. The maximum heat transfer would be obtained if one of the fluids 

were to undergo a temperature change equal to (T, . -T ). the 

h > in c,in * 

maximum temperature difference in the counterflow exchanger. Since the 

heat lost by the hot gas must be gained by the cooler air (neglecting 

losses), as may be seen in Equations 4.9 and 4.10, the fluid undergoing 

this maximum temperature change must have the lower of the two capacity 

rate values, C , . 

min 

q = C . (T. . - T . ) 4.11 

max min h,in c,in 

The effectiveness may now be determined as 

£ - actual heat transfer 
maximum heat transfer 
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€ = CnlTV, J LTV. ~ ^ Vx 

C«v\.n. (JTV^ la - tl ,lrv ) 



CcCq TLvn) 

CL-^aln. (JTV .v.a-'TV yO 



4.12 



From the previous equations, the effectiveness may be shown to be 
independent of the fluid temperatures [2]. 



£ - 




-Cl' 
JL 

a 



C^rOLiry v \3 A. 
C<Aax. 



r v _ C?o>v.ry \ X*J) 

SL V c ^ v - Orv's-n. 



4.13 



where C . and C are, respectively, the smaller and the larger of 
min max y & 

the two capacity rate magnitudes. 

Finally, defining the Number of Transfer Units, NTU (a dimension- 
less expression for the heat transfer size of the recuperator), as the 

ratio of (U A) to C . , where A refers to the total surface area on 

min 

which U is based, the effectiveness may be written as 



£ - 




-U-c^Ont\j 
^ 



C.TVAAC I 

C.cao.'*. 




C— 

0<\Cv< 



) m'tu 



4.14 



A comparison of this equation with effectiveness relations for other 
geometries will lead to the conclusion that the counterflow arrangement 
is indeed useful for gas flow recuperators (for which C . /C 1) 

when an effectiveness in the order of 807 o is to be obtained. 
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B. FLOW - FRICTION THEORY 

In the design of a gas-to-gas heat exchanger, the pressure drop 
characteristics assume an importance equal to that of the heat trans- 
fer characteristics [8]. The designer is generally interested in the 
overall core pressure drop, including entrance and exit effects, rather 
than the core friction pressure drop alone* 




Core Friction and Flow Acceleration 



The core friction pressure loss term, due to viscous shear 

p 



and pressure force, is expressed in terms of the dynamic head 



*9 



the appropriate friction factor, , and the ratio of the length to 



hydraulic radius of one side, — 



A ’ 
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aB 



^ (xe) 



4. 16 



where G is based on the minimum free flow frontal area ( G-p.uw>-- piU_ 2 '- ^ ) 
and p>nv is the average core fluid density ( ). 

Since density changes occur in the core as the fluid temperature 
changes, a term accounting for flow acceleration must be included in 
the pressure drop analysis* The change in momentum due to this fluid 
acceleration may be expressed as 



^ clI ion 



Ac 
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2. Entrance and Exit Effects 
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The entrance pressure drop, due to the free flow area change 
as well as to the irreversible free expansion following the abrupt in- 
let contraction may be written as [8]. 



Px - Bo. - ( p.u.* Ci~ <J~ Z ) -*• p> Kc. 2.0^ 
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Pi - Pa. - 20jp, 1 



4. Kc_ ) 



4.21 



33 



where a constant density is assumed to exist from point "1" to point 
"a 1 . , the contraction, or entrance, coefficient may be calculated 

by methods explained in Ref. [8]. 

Similarly, the exit pressure rise is due to the change in flow 
area (which is identical to the corresponding entrance pressure drop 
term) and the irreversible free expansion and momentum changes follow- 
ing the abrupt exit expansion. This pressure rise is written as 



P\o - 
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Pvs - P2 
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where , the exit coefficient, is evaluated in Ref. [8]. 



3 . Overall Pressure Drop 

The combination of all pressure drop effects results in the 

equation 

G 2 - ' 



Px - P*. = p i 



2.P* / P&. N P* r A 0\ i i \ 
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As it has been assumed that and - p>z the relation for 

the overall exchanger pressure drop in terms of specific volume, where 

may be written as 



Vm. = 

a 



c- 2 vi 

Pi - P 2 r 2 £ 



A Vm Vi 



\ vg. \ r ™ / o \ 

( Wc + 1 - <Z~J v 2C vi • i ^ \ "Ac ~i ' ViClL-cr -Pe.) 
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c. PARAMETRIC APPROXIMATIONS 

In the preceding heat transfer and flow friction analysis, a number 
of important engineering parameters and fluid properties are used in 
fundamental relationships. Usually presented only in tabular or graphi- 
cal form, these parameters may not be directly obtained in a computer 
analysis. To overcome this problem, functional relationships approxi- 
mating the tabulated values must be made a part of the computer program. 

For the program [1], such relationships were formulated by making a 
"least squares fit" of the curves representing tabulated data of the 
required parameters. The resulting functional relationships, of the 
form 

- A; + A 2 X- t AjK. 2 +■ v Av-X}' 1 4.26 

where: P(x) = desired parameter 

x = variable on which the 
parameter is based 

A^ = polynomial coefficients 
resulting from least 
squares fit 

i-1 = order of resulting 
polynomial 

are used extensively in the computer analysis. Table II specifies the 
parameters for which such polynomial approximations are used, the 
variable on which each is based, and the order of the polynomial. 

Table II indicates that two parameters for which functional relation- 
ships must be obtained are dependent on the Reynolds number , (4 /jjl. 

The Fanning friction factor, § , is used in equation 4.15 and subsequent 

relations to obtain the pressure drop due to core friction through the 
recuperator. The Colburn- j factor, j, has been obtained experimentally 
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for a number of recuperator geometries [8] and is used to obtain the 
convection heat transfer coefficient, h, as follows: 

This heat transfer parameter is defined as 

2/3 

j - St Pr ' 4.27 

Since j and Pr are determined from functional relationships, the 
Stanton number may be calculated as 

-2/3 

St = j Pr 4.28 

Thus, the convection heat transfer coefficient, used in equation 4.5 
and subsequent relations, is obtained as 

h = St 0 ^ p u 4.29 

Experimental investigations of ^ and j provide curves similar 
to those presented in Figure 7. As indicated, the laminar flow regime 
(Re ^ 1000) produces linear relationships for the ideal equilateral 

triangular passage such that 



f = 13.33/Re 


4.30 


j = 2.64/Re 


4.31 



The turbulent region (Re > 5000) data may be closely approximated by 

the linear relationship resulting from circular tube flow as 

f = 0.046 Re'°” 2 4.32 

j = 0.023 Re' 0 * 2 4.33 
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Re 

FIGURE 7 

EXPERIMENTAL CURVES FOR -f AND j 

The transition region (1000 i Re £. 5000) is also based on experi- 
mental data and must be approximated by a polynomial equation for 
computer analysis. The polynomial coefficients for the resulting 
equations approximating this transition region, as well as those for 
the other parameters in Table II, are presented in Appendix B. 
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TABLE II 



FUNCTIONAL RELATIONSHIPS 





Parameter 


Function Of 


Order Of 
Pol ynomi a 1 




Fanning friction 
factor 


Reynolds number 


10 


j 


Colburn-j factor 


Reynolds number 


10 


Pr 


Prandtl number 


temperature 


9 


JU. 


vl scos I ty 


temperature 


9 


C P 


specific heat 


temperature 


5 


1 


entha 1 py 


temperature 


6 


P re 1 


relative pressure 


temperature 


8 


T 


temperature 


relative pressure 


8 


T 


temperature 


enthalpy 


9 
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V. RESULTS 



A. OBJECTIVES 

For many years engineers have designed heat exchangers by first 
assuming core dimensions and then calculating mass flow rates, heat 
transfer rates, and pressure drops which would occur if fluids were 
flowing through the assumed core. Extensive investigations of various 
heat transfer surfaces offer the designer enough latitude in surface 
geometry to successfully meet the heat transfer and pumping power re- 
quirements for a specific application. 

Texts frequently refer to the possibility of optimizing the inter- 
nal core dimensions of a heat exchanger once a particular surface is 
specified, yet very little research has been published in this area. 
Such a technique forms the basis for this study. 

The computer program [1] used in the study simulates the open gas 
turbine regenerated cycle under constant load conditions. All adjust- 
ments in internal exchanger geometry are reflected directly in terms of 
the cycle thermal efficiency and the overall core size necessary to 
obtain that efficiency. As is outlined in Appendix C, any computer 
simulation of the gas turbine cycle requires set tolerance limits on 
various cycle design points which must be met. Early investigations of 
general recuperator trends in the study involved a tolerance limit of 
3.3% for the desired thermal efficiency of 0.34 and a limit of 0.2% for 
the desired horsepower output of 600. More detailed analysis (as pre- 
sented in Figures 15 through 22 and Tables V through VII) required 
tighter tolerances of 1.45% for the thermal efficiency and 0.097 o for 
the horsepower. The use of two tolerance levels slightly affects out- 
put data for the same core, as may be noticed. 
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B. GENERAL RECUPERATOR CHARACTERISTICS 

Trends in external dimensions of the plate-fin recuperator due to 
changes in internal spacings may be noted by investigating a number of 
different cores. A sufficient number of cores may be obtained by pair- 
ing the plate spacings presented in Table I. For convenience a parti- 
cular core will be referred to in terms of the ordered pair (b ,b, ) where 

c h 

b represents the cold side spacing and b. represents the hot side spac- 
c h 

ing of the core . 

Due to addition of fuel in the combustor, the mass flow rate through 
the recuperator hot side is larger than the cold side flow rate by a 
factor of the fuel to air ratio, r £/ a * The hot side mass flow rate must 
be accommodated, then, by an increase in fluid density, flow area, or 
fluid velocity (m - p A v) . Since fluid density decreases as tempera- 
ture increases the engineer is left with only two alternatives. 

A velocity increase will not only accommodate the increased mass 

flow rate, but it will also increase the heat transfer rate per unit of 

■ » 

surface area by a factor something less than the first power of the 
velocity. The friction power expenditure, however, will increase by a 
factor between the square and the cube of the velocity [8]. It is this 
effect which serves as the basis for the design restriction that the hot 
side plate spacing never be less than the cold side spacing. 

General recuperator volume trends may be noted in Figure 8 for five 
specific cores. Each point on every curve represents the external core 
dimensions which will satisfy the specified design points. As might be 
expected, from the definition of the compact heat exchanger, larger cores 
are required as one moves from the more compact to the larger spacings; 
for any one frontal area, more length, and thus more heat transfer area, 
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is required as spacings increase. Excluding the extremely compact core, 
(0.1, 0.1), points of minimum volume for each core fall below a frontal 
area of five square feet. Furthermore, every core requires a frontal 
area of at least two square feet before the specified cycle design points 
are met. 

A more detailed view is presented in Table III, in which the minimum 
volume of the (0. 1,0.1) core is seen to fall within a frontal area range 
between 30 and 40 square feet. Definitely, the minimum volume of such a 
core does not present the optimum package. Similar characteristics were 
noted only for two other extremely compact cores, (0.1, 0.125) and (0.125, 
0.125). 

An increase in length over a region of increasing frontal area is 
noticeable in all cores except those having a cold side spacing of 0.1 
inches. Occurring at a frontal area greater than that for the minimum 
volume core, this effect becomes more apparent as core spacings increase. 
Table III indicates that the mass flow rates for these cores decrease as 
frontal area increases. Not shown is the fact that the recuperator hot 
and cold side Reynolds numbers also decrease with the mass flow rate. 

The first knee of each curve (that point at which L begins to increase 
with A^ ) occurs when the two Reynolds numbers approach a magnitude of 
2000. Reference to Figure 7 indicates a decrease in the Colburn-j fac- 
tor with decreasing Re at this point. Evidently, the resulting decrease 
in the convection heat transfer coefficient, h (which is proportional to 
j), is accompanied by an increase in heat transfer area (as L increases 
with A^) in order that the desired cycle thermal efficiency be main- 
tained. This effect continues as j decreases until the second knee at 
which the two Reynolds numbers have approached a magnitude near 1000. 
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At this point, since j begins to increase with increasing frontal area, 
less heat transfer area is required, and L begins to decrease. 

The effect is not noticeable for those cores of low plate spacings, 
as both the hot and cold side Reynolds numbers of such cores fall in the 
laminar region for all frontal areas . Reynolds numbers for cores of low 
cold side spacing and large hot side spacing (i.e., (0.1, 0,5)) however, 
tend to have a maximum Reynolds number near 2000 at low frontal areas, 
followed by a rapid drop as frontal area increases. The slight length 
increase required by such cores is not noticeable due to the magnitude 
of the length increments used at this point of the investigation. 

Relative pressure drops on both sides of all cores decrease to zero 
as frontal areas increase. In the region of minimum volume (A^ 4 - 5 sq 

ft), the cold side relative pressure drop appears to be negligible for 
all cores except that which combines a small cold side spacing and a 
large hot side spacing. The hot side relative pressure drop, on the 
other hand, appears to be high only for those cores in which the two 
plate spacings are of equal magnitude. 

C. MINIMUM VOLUME AND WEIGHT CORES 

Figure 9 indicates dimensional requirements for various cores in the 
minimum volume region. As is to be expected, required lengths, and there- 
fore volumes, increase as core spacings increase. The effect of forming 
cores of unequal spacings as compared to equally spaced cores is presented 
in Figure 10. 

To this point a trend has developed in which smaller spacings re- 
quire smaller cores. Based on this trend, one might incorrectly reason 
that the minimum volume for a given hot side spacing may be obtained 
simply by decreasing the cold side spacing to a minimum allowable magni- 
tude. In fact, for the cores of 0.25 inch hot side spacing this appears 
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to be the case as b is reduced from 0.25 inches to 0.1 inches. How- 

c 

ever, the trend presented by those cores for which b^ is 0.5 inches as 
b^ decreases should be noted. A decrease from a cold side spacing of 
0.5 inches to 0.25 inches reduces core volume; a further decrease to a 
b^ of 0.1 inches produces a core of extremely large volume. This core, 
which combines the large length and frontal area requirements of the two 
range limiting spacings, results in the largest minimum volume point re- 
corded of all cores investigated. 

Similar analysis shows the effect, presented in Figures 11 through 
13, of increasing the hot side spacing for a constant b^. For all cores 
presented, the b^ increase results in increased frontal area and length 
requirements; thus, the minimum volume for each series of constant cold 
side spacing cores increases as the hot side spacing is increased. 

Recuperator dimensional and performance characteristics at minimum 
volume points are presented in Table IV ; each set of cores listed con- 
sists of a constant b^ with increasing b^. The optimum core with re- 
spect to volume in each set occurs at or near that core configuration 
for which b^ is twice b^. The optimum core with respect to weight for 
each set consists of the same, or a slightly higher, b^ as the optimum 
volume core. Such a trend might lead one to believe a rule of thumb 
often used by manufacturers that the core hot side spacing be twice the 
cold side spacing. Analysis will show this to be an inadequate design 
procedure. Further investigation in this area, however, should be pre- 
ceded by mention of the performance characteristics of the minimum core 
points as tabulated in Table IV. 

For any one set of cores, the cold side relative pressure drop in- 
creases with increased b, , while the hot side relative pressure drop 

h 
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decreases. This effect is presented graphically for one set of cores in 

Figure 14. The large magnitude of the hot side relative pressure drop 

for cores of b, t b is to be expected since the hot side mass flow rate 
h c 

must be accommodated by an increase in either flow area or velocity over 

that of the cold side flow. Since there is little or no area increase 

(b^ X b^), velocities on the hot side of such cores must be higher, 

yielding relatively large pressure drops. As mass flow rates generally 

fall and flow areas increase as b, is increased for each set of cores, 

h 

the hot side relative pressure drop falls „ Evidently, the cold side 

relative pressure drop is affected more by the general length increase, 

as b, increases for each set, than by the decreasing trend in the mass 
h 

flow rate. Thus, the cold side relative pressure drop rises with b^ in 
each set of cores. 

D. OPTIMUM CORES 

Concentrating once again on dimensional aspects. Figures 15 through 
19 represent the effect of increasing b^ for a constant b^ in the region 
of optimum core spacings with respect to volume and weight. It is ap- 
parent from these plots that an optimum design procedure is not that for 

which b, “ 2b . Tables V and VI present the characteristics of cores 
h c 

optimized with respect to volume and weight; Table VII, which consists 

of performance data for cores of b, = 2b , is presented for comparison. 

h c 

A comparison of Tables V and VI indicates that the minimum volume core 

for any one b involves a higher mass flow rate and less length than the 
c 

corresponding minimum weight core. This results in a lower cold side 
and higher hot side relative pressure drop through the minimum volume 
core as compared to the minimum weight core for any set b^ . 
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Further investigation of the Tables indicates that the minimum weight 

core for any requires a larger b^ than that for the minimum volume 

core. This effect (noticed earlier in Table IV) is to be expected, as 

the minimum weight core represents a compromise between a low core volume 

and a low overall core density (which decreases as spacings increase). 

Differences between the two sets of optimized cores and those cores for 

which b, = 2b are illustrated in Figure 20. 
h c 

Allowing a tolerance limit of 0.004 inches either side of the optimum 
hot side spacing for any given cold side spacing as presented in Tables 
V and VI, smooth optimum spacing curves may be drawn, as shown. In the 
range of cold side spacings between 0.1 and 0.225 inches a linear re- 
lationship closely approximates each of the two optimum spacing curves, 
as 

b u = 0.86 b + 0.12 (0.1 1 b 1 0.225) 5.1 

he c 

for optimum volume spacings and 

b - 1.16 b + 0.12 (0.11b £ 0.225) 5.2 

he c 

for optimum weight spacings. 

Outside this range of cold side spacings the optimum core plots 
tend to level off, thus tending to keep the cores within the limits of 
compactness of the compact heat exchanger region. 

Since the engineer is usually at liberty to choose both plate spac- 
ings, he must have an idea of the relative volumes and weights between 
the optimum cores which he may use. Figure 21 presents volume character- 
istics of optimized cores from Table V with respect to the minimum volume 
core (0.125, 0.225). As may be seen, if the engineer expects to use a 
cold side spacing equal to or greater than 0.175 inches, he must accept 
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the fact that his optimum core will have a volume at least 207* greater 

than the minimum available* Similarly,, Figure 22 represents the weight 

characteristics of cores from Table VI with respect to the minimum weight 

core of (0.175, 0*332). Although the variance in weight is not as great 

as that for volume, on a percentage basis 9 it may be seen that the enginee 

must expect weights at least 57* greater than the minimum available if he 

wishes to use cold side spaeings other than those between 0.125 and 0.235 

inches. The very large increase in relative weight as extremely compact 

cores (b^ ^ 0.1 inches) are optimized should also be noted. 

For the design points of this study, if the engineer were free to 

optimize any non-dimensional parameter (as mass flow rate or relative pres 

sure drops) while keeping within 107* of the minimum volume possible, refer 

ence to Figure 21 and the use of equation 5-1 would indicate a freedom of 

b choice between 0.1 and 0.15 inches and b. choice between 0.206 and 0.28 
c h 

inches. A similar design procedure with a limit of 57* over the minimum 

weight core (using Figure 22 and equation 5-2) would yield a b range be- 

c 

tween 0.13 and 0,225 inches and a b, range between 0.27 and 0,38 inches. 

h 

A reference back to Table IV will indicate that the minimum weight and 
minimum volume cores for all cores investigated do fall within the indi- 
cated ranges of spaeings. If both volume and weight of the core are 
limited within the specified 107* and 57* 9 the range of allowable cold side 
spaeings is further reduced to between 0.13 and 0.15 inches, and the allow 
able hot side range to 0.27 to 0.28 inches. 

Although the design limitations outlined are applicable only to the 
specific core geometry studied, the procedure may be applied to all simi- 
lar geometries. As indicated, when core and weight tolerances are speci- 
fied, the engineer does have some freedom in choosing internal core 
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dimensions based on other requirements. In all applications for which 
the allowable volume or weight of the gas turbine is at a premium, it 
seems that the engineer would be wise to follow a design procedure for 
the heat exchanger, as outlined, which optimizes the limited quantity, 
before moving on to further core optimization with respect to perform- 
ance. 
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FIGURE 8 GENERAL EFFECT OF PLATE SPACINGS ON CORE DIMENSIONS 



GENERAL RECUPERATOR PERFORMANCE AND DIMENSIONAL TRENDS 
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FIGURE 9 CORE DIMENSIONS FOR 
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FIGURE 10 EFFECT OF UNEQUAL SPACINGS ON CORE DIMENSIONS 




RE DIMENSIONS OF INCREASING b. FOR b„ = 0.1 INCHES 
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FIGURE 12 EFFECT ON CORE DIMENSIONS OF INCREASING b. FOR b - 0.2 INCHES 
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FIGURE 13 CORE DIMENSIONAL BEHAVIOR AT LARGE SPACINGS 




THE EFFECT ON PERFORMANCE AND DIMENSIONAL CHARACTERISTICS 
OF INCREASING b L FOR CONSTANT b 
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TABLE IV 
( contl nued ) 
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11.96 373.5 .018 .035 .743 21467 

15.44 415.0 .025 .025 .740 21381 



TABLE IV 
(continued) 
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FIGURE 14 RELATIVE PRESSURE DROP AS A FUNCTION OF FOR b * 0*1 INCHES 
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L (ft) 

FIGURE 15 EFFECT OF INCREASING ON CORE DIMENSIONS IN OPTIMUM REGION 

FOR b » 0.1 INCHES 
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FIGURE 16 EFFECT OF INCREASING b h ON CORE DIMENSIONS IN OPTIMUM REGION 

FOR b r = 0.125 INCHES 
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L (ft) 

FIGURE 17 EFFECT OF INCREASING b h ON CORE DIMENSIONS IN OPTIMUM REGION 

FOR tv a 0.15 INCHES 
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L (ft) 

FIGURE 18 EFFECT OF INCREASING b h ON CORE DIMENSIONS IN OPTIMUM REGION 

FOR b_ - 0.175 INCHES 
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FIGURE 19 EFFECT OF INCREASING b h ON CORE DIMENSIONS IN OPTIMUM REGION 

FOR b - 0.2 INCHES 



CHARACTERISTICS OF MINIMUM VOLUME CORES 
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CHARACTERISTICS OF MINIMUM WEIGHT CORES 
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FIGURE 20 OPTIMUM RECUPERATOR CORE SPACINGS WITH RESPECT TO WEIGHT AND VOLUME 
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FIGURE 21 RELATIVE VOLUMES OF OPTIMIZED CORES WITH RESPECT TO THE MINIMUM VOLUME CORE 
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FIGURE 22 RELATIVE WEIGHTS OF OPTIMIZED CORES WITH RESPECT TO THE MINIMUM WEIGHT CORE 



APPENDIX A 



THE REGENERATED BRAYTON CYCLE 



A. THERMODYNAMICS 

Ideally, the basic open cycle consists of the four processes of 
isentropic compression, constant pressure heat addition, isentropic 
expansion, and constant pressure heat rejection [4]. The ideal regener- 
ated cycle differs only in that the heat addition process is accomplished 
through the use of a recuperator as well as the standard combustor. With 
these points in mind, the ideal cycle (refer to Figure 23) may be des- 
cribed using the equalities: 



P 2a 


P 3a P 4a 


A.l 


P 5a 


P 6a 


A. 2 


T 2a 


= T, 
6a 


A. 3 


T 3a 


= T 

5a 


A. 4 



The most obvious differences between the actual and ideal cycles^ 
from the thermodynamic viewpoint, are the increases of entropy and drops 
in pressure due to actual working conditions. These differences may be 
best explained by taking the cycle in its respective stages. 

1 . Compression 

State point n 2 n lies to the right, at an increased entropy, of 
ideal point 9C 2a n due to the irreversibility of the compression process; 
both points, however, lie on a line of constant pressure, since the 
compressor outlet pressure is specified by working conditions. The 
resulting equality is 
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("a" Indicates Ideal 
state points) 




FIGURE 23 THE ACTUAL AND IDEAL REGENERATED BRAYTON CYCLES 
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2 . 



Heat Addition 



P 



2a 




A. 5 



Two changes are apparent when one compares cycle conditions at 
the recuperator exit point* First, the actual pressure level is below 
that of the ideal, due to the pressure drop as the air flows through 
the recuperator* 



P~ < P~ 

3 3a 



A. 6 



Second, the temperature of the actual case is below that for the ideal, 
since less heat is transferred than for the reversible process. 



T < T 
3 3a 



A. 7 



In order to raise the temperature further, the air flows through a 
combustor. The inlet temperature of the turbine is specified by its 
material properties. 



T' - T 
4a" X 4 



A. 8 



However, point ,f 4 ff does lie at an entropy state greater than point "4a n 
due to the irreversibilities involved between the compressor and the 
turbine. Furthermore, an additional pressure drop occurs as the air 
flows through the combustor. 



P, ^ P, 

4 4a 



A. 9 



3 . Expansion 

An additional entropy increase results from the expansion 
process in the turbine, moving point n 5 n to the right of the ideal state 
point. As the turbine outlet pressure is not specified, and a further 
pressure drop is expected as the gases pass through the recuperator, 
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the pressure at point M 5 H is necessarily greater than that for the 
ideal cycle. 




A. 10 



4. Heat Rejection 

As indicated, the pressure drop experienced by the hot gases 
in passing through the recuperator results in the pressure at point ,, 6 ,, 
being equal to the ideal pressure, since this point is considered to be 
the exhaust to atmospheric conditions. 



Finally, since less heat is transferred from the hot gases to the in- 
coming cold air as is ideally possible, the temperature at point 1 *6 ,f 
is greater than that for the ideal case. 



B, CYCLE SYSTEMS 

An understanding of the efficiencies of each system is necessary 
if one is to work with the Brayton cycle. This may best be accomplished 
by looking at each system process individually. It must be noted that 
the assumption of steady flow and no potential energy change of posi- 
tion is made for the cycle. More detailed information may be obtained 
in Refs. [3] and [5] . 

1 . Compressor 



P 6a = P 6 = 14.696 psi 



A. 11 




A. 12 



= compressor efficiency 



Qc 



work of isentropic compression = X la 1 2a 
actual internal work -i 



A. 13 
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2 . 



Recuperator 

C = recuperator effectiveness 



€. = 



actual heat transfer 

maximum possible heat transfer 



£ ~ - ^,out) __ Crn <- ^Tc .ovjl -'"Pc.in) * ^ 

C ^C^v)tCLV.rv C% i v.tx -T, c,V a') K ' ^ C^Vi ,lr\ -~^caa) 



A more detailed explanation of recuperator effectiveness is 
presented in Section IV 0 

3. Combustor 



combustor efficiency 

actual energy received by gases 

energy released during complete combustion 

U^JT.^olMa - is - JjZSfe. ii A>15 

This definition may be more easily understood when an energy bal- 
ance based on one pound of air and £? pounds of fuel per pound of 

air entering the combustor is made. 

* r f/a. 4.X + s + A ‘ 16 



^cor-r^ = 
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l ~ r ;/cc 



^1 4 r ^/(X ^a) 



of the three entering energy terms, it can readily be seen that i^ (the 
enthalpy of the incoming air) and r f/ a *-f (the fuel to air ratio times 
the enthalpy of the fuel) represent molecular internal energy and flow 
work. The third term represents the chemical energy stored in the fuel, 
r f/a^l* 5 s i nce complete condensation of water vapor in the air flow will 
not occur, the lower heating value of the fuel (which disregards the 
latent heat of condensation), q^, is usually used in gas turbine anal- 
ysis [3] . 

The fuel to air ratio, a very important parameter in the gas tur- 
bine analysis, may now be shown to be 



4. Turbine 



at 

at 



= turbine efficiency 

actual internal work 

S ' — 1 ■■ 

work of isentropic expansion 
_ i-4 - 

i - i^ou 



A. 17 



A, 18 
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5. 



Additional Calculations 



ysxs , 



Three final calculations are frequently useful in a complete anal- 

Overall Thermal Efficiency 

- overall thermal efficiency 



O-U-L ■ 



net work 
heat supplied 



o . u cl ~ 



A. 19 



Horsepower 

Using the conversion factor 2544 BTU = one horsepower- 
hour, the equation for the horsepower output may be seen to be 



hp = (Si ;/oJ C q 1 i ) a . 20 

2544 

Specific Fuel Consumption 



SFC = 


fuel burned 




net work 




SFC = 


Crf/cO 

HP 


A. 21 


SFC = 


2544 


A. 22 


*>1 
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APPENDIX B 



COMPUTER PROGRAM ESSENTIAL FEATURES 
A. GENERAL ASPECTS 

The computer program referred to in Ref. [1] has been greatly ex- 
panded by Professor P. F. Pucci since its original use, and may now be 
used to investigate parallel and crossflow, as well as counterflow, 
recuperator characteristics* This counterflow analysis is based on 
design parameters for the Orenda OT-4 gas turbine only; the program, 
however, is applicable to any gas turbine recuperator analysis, as the 
desired design points are introduced as input data. A basic under- 
standing of the program logic may be obtained by referring to the simpli- 
fied flow chart, Figure 25; a sample program output is presented in 
Figure 28. 

The program is written on the basis that changes in recuperator core 
geometry (plate spacings and fin configurations) and dimensions (frontal 
area and length) are reflected directly in the overall cycle thermal 
efficiency, rather than in the often used design variables of heat ex- 
changer pressure drop and heat transfer characteristics. Although the 
program may be used to obtain data for a number of conditions (such as 
the computation of: (1) the maximum cycle thermal efficiency given a 
recuperator of specified volume, (2) recuperator dimensions necessary to 
obtain a specified thermal efficiency and, (3) general outputs for a 
simple open cycle, to mention only a few) only two output types were used 
extensively in the analysis. 
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The first, Output Type III, computes the required recuperator 
length for a specified frontal area and internal geometry; furthermore, 
the horsepower and cycle thermal efficiency are specified, within pre- 
determined tolerance levels* This output is extremely useful in provid- 



range of frontal areas* Information on the effect of changes in internal 
core dimensions, may be obtained with the use of Output Type VIII. This 
routine, similar to Output Type III, computes the required length of a 
recuperator over a range of frontal areas* In this procedure, however, 
the cold side plate spacing remains constant, as before, whereas the hot 
side plate spacing is incremented over a desired range. Flow charts for 
these routines are presented in Figures 26 and 27. 

B . PARAMETRIC APPROXIMATIONS 

As explained in Section IV, the computer program makes use of poly- 
nomial approximations to obtain a number of fundamental parameters. The 
polynomial coefficients used in this analysis will be presented, should 
a need for them arise in the future; unless otherwise indicated, the 
polynomial is of the form 



ing an overall picture of the desired recuperator's behavior over a wide 




+ A . x 



1 



i -1 



B.l 



where: P(x) = desired parameter 



x = 



variable on which the parameter is based 



A^ = polynomial coefficients 
i-1 = order of polynomial 
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Fanning Friction Factor (function of Re; 1000 Re ^ 5000) 



A x = -1.311518618929068 x 
A 2 = 7.522772995845896 x 
A 3 = -1.550552244236629 x 
A, - 1.730155685156651 x 
A 5 = -1.189041200805969 x 
A, = 5.329916594794135 x 

D 



10' 1 A ? = -1.590750036125343 x 

10" 4 A = 3.138452606608756 x 

O 

10 -6 A 9 = -3.933692644883426 x 

10' 9 A 1q = 2.838094432805311 x 

10' 12 A 1X = -8.976537263266828 x 



10 

10 

10 

10 

10 



-19 

-23 

-27 

-31 

-36 



Colburn- j Factor (function of 

A 1 = 1.644630125351775 x 

A 2 = -6.796500408006909 x 

A 3 = 1.238993655152145 x 

A, = -1.292429783097030 x 
4 

A 5 = 8.543160388163581 x 

A, = -3.740734237674230 x 
6 



Re 

10 

10 

10 

10 

10 

10 



; 1000 * Re 5000) 

_1 A ? = 1.100442463486636 x 

A. = -2.152007282686803 x 

O 

A 9 = 2.683333339022575 x 

A = -1.930653368061576 x 

" 13 A = 6.099637443396928 x 

-16 



10 

10 

10 

10 

10 



-19 

-23 

-27 

-31 

-36 



Prandtl Number (function of temperature) 



A = 0.7865630968 

A 2 =-0.2871584338 x 10‘ 

A 3 = 0.3666074673 x 10 

A. =-0.1373281188 x 10 ‘ 
4 

A 5 = 0.1569904076 x 10 



-€> 



-11 



A £ = -0.1973361281 x 10 

D 

A ? = 0.1041204710 x 10 

A 0 = -0.2047252924 x 10 

O 

Ag = -0.1538501631 x 10 
A = 0.7960530427 x 10 



■14 



-17 



■21 



■25 



-29 



Viscosity (function of temperature) 



- 10 , 



JU. 
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(as a function of temperature) of coefficients: 



80 



A = -0.4132588012 x 10 

A 2 = 0.2706898585 x 10' 

A 3 = -0.1318733093 x 10 

A, = 0.7538864121 x 10 

4 

A 5 = 0.7221455415 x 10 



-5 



-12 



■12 



A 6 = 



A 8 

A 9 " 
A 10= 



-0.5752187450 x 10 



0.2117993215 x 10 



-15 



■18 



-0.3862860359 x 10 



0.2782105119 x 10 



0.6763288488 x 10 



-22 



-26 



-32 



Specific Heat (function of temperature) 

,0 



A = 0.24913875 x 10 

A 2 = -0.48196782 x 10 

A = 0.68100588 x 10‘ 

3 



-4 



A 4 = 
A 5 " 
A 6 " 



-0.21358038 x 10 



0.48224330 x 10 



-10 



-15 



0.47985777 x 10 



■18 



Enthalpy (function of temperature) 

2 



A = -0.97746770 x 10 
A 2 = 0.25257199 x 10 
A 3 = -0.31167200 x 10‘ 



0 



A, = 0.30219636 x 10 
4 

A 5 = -0.94272843 x 10‘ 
A 6 = 0.11868992 x 10' 
A ? = -0.32615268 x 10' 



■11 



■14 



■19 



Relative Pressure (function of temperature) 
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Temperature (function of enthalpy) 



The temperature is computed in the standard polynomial form, as a 
function of x, where 

x = log(10.0 * temperature) B.4 
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FIGURE 25 PROGRAM FLOW CHART 
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FIGURE 27 OUTPUT VIII FLOW CHART 
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UTPUT TYPE THREE 
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APPENDIX C 



ERROR ANALYSIS 



It is a basic fact that all engineering investigations are subject 
to inaccuracies. As in any other study, a computer analysis is in- 
herently in error due to: (1) the basic assumptions involved to facili- 

tate the problem solution (as neglecting plate heat transfer resistance 
in equation 4.7) and (2) the approximations made in order to use pre- 
viously compiled engineering data (as the polynomial curve fittings dis- 
cussed in Appendix B) . A complete examination of all the errors 
introduced by these factors would be not only time consuming, but quite 
pointless. 

The most significant error in this computer analysis arises from 
the use of tolerance limits on necessary output parameters and finite 
increments on input parameters. This investigation necessitates the 
use of two sets of tolerance levels. The first. Set I, involves a 
tolerance limit of 3.3% (0.3389 ^ l 0.3411) on the desired 

thermal efficiency of 0.34 and a limit of 0.27> (599.0 £ HP - 601.0) 
on the desired 600 horsepower. Data for these limits is calculated at 
frontal area increments for the recuperator core of 0.02 sq ft and 
length increments of 0.01 ft. Further investigation, in the area of 
optimization, requires tighter tolerances of 1.45% (0.3395 — - 0.3405) 

on the thermal efficiency and 0.097, (599.5 - HP - 600.5) on the 

horsepower. Data for these levels, designated as Set II, is supplied 
at frontal area increments of 0 o 005 sq ft and length increments of 
0.005 ft. Data from Set I is reproduced in Figures 8 through 14 and 
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Tables III through IV; Set II data is presented in Figures 15 through 
22 and Tables V through VIII. 

As there are no performance characteristics available for the 
majority of recuperator cores investigated in this study, no judgement 
can be made as to the errors introduced by the computer simulation. An 
investigation of variations in core parameters as presented in Table VIII, 
however, yields an idea as to the relative error introduced by the neces- 
sary use of finite tolerances and increments. 

It should first be noted in Table VIII that three data sets are 
presented for each core. Set III involves the use of tolerance levels 
as specified for Set II, with reduced frontal area and length increment 
magnitudes of 0.002. Variations between Sets II and III are so small 
that the added computer calculation time introduced with the use of Set 
III over-rides the slight benefits due to the smaller increments. 

Variations in all outputs except dimensional and weight parameters 
are of such an order that the error in computing these quantities seems 
negligible. The most noticeable differences between output sets arise 
in the frontal area and length requirements and the subsequent weight 
and volume quantities. The effect of tolerance levels on dimensional 
requirements is more noticeable in Figures 29 and 30. 

As may be seen, for any given frontal area, the tighter tolerance 
of Set II produces a core of greater length than Set I. A comparison 
of the two plots results in the fact that this effect is more notice- 
able at larger core spacings. For all cores tested, length variations 
between sets decrease as frontal areas increase and less core length is 
required. It is difficult to establish a trend as to the frontal areas 
of minimum volume cores for the two sets. Although Table VIII indicates 
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that Set II produces minimum volume cores of greater frontal area than 
Set I, this effect does not hold for all cores investigated. However, 
all Set II minimum volume cores of less frontal area than the corres- 
ponding Set I minimum volume core require a much greater length than 
that for the Set I core. Thus, the overall volume and weight trends 
between cores for each specific set seem to hold. 

The most important questions concerning relative error of the analy- 
sis have been answered. There is evidently little error in the perform- 
ance characteristics output. Dimensional comparisons between the two 
sets result in a maximum difference in length of 0.1 ft for any given 
frontal area and 0.1 sq ft in area for any minimum volume core. Most 
important, due to similar trends in the two sets, the data presented 
may be assumed to provide specific patterns in plate-fin counterflow 
recuperator core dimensions and weight requirements, while the actual 
magnitudes presented may be assumed to be within a reasonable range of 
those which might actually be required. 
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THE EFFECT OF TOLERANCE LIMITS ON MINIMUM VOLUME DIMENSIONS 

OF SPECIFIC CORES 
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91 



FIGURE 29 EFFECT OF TOLERANCE LIMITS ON (.125,. 225) CORE 
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FIGURE 30 EFFECT OF TOLERANCE LIMITS ON (.4, .5) CORE 




BIBLIOGRAPHY 



1. Carpenter, So P„ Computer Program for Sizing Gas Turbine Counter - 
flow Recuperators , Master's Thesis, Naval Postgraduate School, 
Monterey, 1967 o 

2. Chapman, A. J. , Heat Transfer , 2d ed., The Macmillan Company, 
1967. 

3. Dusinberre, G. M, and Lester, J. C., Gas Turbine Power , 2d ed. , 
International Textbook Company, 1958. 

4. Faires, V. M. Thermodynamics , The Macmillan Company, 1962.. 

5. Howard, C. P M "'Thermodynamics and Fundamentals of the Gas Turbine 
Cycle", and "Cycle Calculations and Variations", Gas Turbine 
Engineering Handbook , Gas Turbine Publications, Inc., 1966. p 1-28. 

6. Howard, C. P., "Heat Exchangers", Gas Turbine Engineering Hand - 
book , Gas Turbine Publications, Inc., 1966. p 183-190. 

7. Hryniszak, W. , Heat Exchangers , Butterworth *s Scientific Publica- 
tions, 1958. 

8. Kays, W. M. and London, A. L. , Compact Heat Exchangers , 2d ed., 
McGraw-Hill Book Company, 1964. 

9. London, A. L. , "Compact Heat Exchangers", Mechanical Engineering , 
v. 88, p 47-51, May 1964, p 31-34, June 1964, p 33-35, July 1964. 

10. Quan, D. , "The Design and Development of the Orenda OT-4 Gas 

Turbine", Trans. ASME, Journal of Engineering for Power , v. 88, 
p 117-126, 1966. 



93 



INITIAL DISTRIBUTION LIST 



No* Copies 

1. Defense Documentation Center 20 

Cameron Station 

Alexandria, Virginia 22314 

2. Library, Code 0212 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Naval Ship Systems Command, Code 2052 1 

Department of the Navy 

Washington, D* C, 20360 

4. Professor P. F* Pucci, Code 59Pc 5 

Mechanical Engineering Department 

Naval Postgraduate School 
Monterey, California 93940 

5. LTJG B. V, Burrow, USN 2 

533 McGilvra Blvd. , East 

Seattle, Washington 98102 

6. Mechanical Engineering Department 2 

Naval Postgraduate School 

Monterey, California 93940 



94 



Unclassified 

Security Classification 



DOCUMENT CONTROL DATA - R & D 

(Security classification of title, body of abstract and indexing annotation must be entered when the overall report Is c lassltlcd) 



l, ORIGINATING activity (Corporate author) 

Naval Postgraduate School 
Monterey, California 93940 



2*. REPORT SECURITY CLASSIFICATION 



Unclassified 



2b. GROUP 



3 REPORT TITLE 



The Effect of Recuperator Geometry on a Regenerated Brayton Cycle 



4. DESCRIPTIVE NOTES (Type of report and.inclusivc dates) 

Master's thesis; June 1969 

5 au THORIS) (First name, middle initial, last name) 

Barry Vaile Burrow 



6 REPORT DATE 



June 1969 



7a. TOTAL NO. OF PAGES 



95 



7b. NO. OF REFS 

10 



8a. CONTRACT OR GRANT NO 



b. P ROJEC T NO 



9a. ORIGINATOR’S REPORT NUMBER(S) 



9b. OTHER REPORT NO(S) (Any other numbers that may be assigned 
this report) 



d. 



10 DISTRIBUTION STATEMENT 

Distribution of this document is unlimited 



1 1 . 



SUPPLEMENT ARY NOTES 



12 SPONSO RING Ml LI T AR Y ACTIVITY 



Naval Postgraduate School 
Monterey, California 93940 



A computer simulation of the Brayton cycle is used to develop a design 
procedure with respect to minimizing volume and weight for the counterflow 
plate-fin recuperator. Based on the Orenda 0T-4 600 horsepower gas turbine, 
recuperator performance and dimensional characteristics are presented for an 
idealized equilateral plate-fin recuperator core matrix. The effects of 
plate spacings varying between 0.1 inches and 0.5 inches on recuperator 
performance characteristics are presented over a wide range of core frontal 
areas. Specific trends toward minimum volume and minimum weight plate-fin 
recuperator cores are discussed in detail. 



DD 



FORM 

1 NOV 85 



1473 



(PAGE 1) 



Unclassified 



S/N 01 01 -807-601 1 



95 



Security Classification 



A- 31400 



Unclassified 



Security Classification 



1 4 

KEY WORDS 


LINK A 


LINK B 


LINK C 


ROLE 


W T 


ROLE 


W T 


ROLE 


W T 


Gas Turbine 

Regenerated Brayton Cycle 
Recuperator 












1 

1 



DD 



FORM 

t NO V 69 



1473 (BACK) 



96 



Unclassified 



/N 0101-307-69?! 



Security Classification 



A- 3 \ 409 




SHELF BINDER 



: Syracuse, N, Y c 

— Stockton, Calif. 

! 


















thesB88518 

° f recupe,ator geometry on re 




3 2768 002 08802 3 
DUDLEY KNOX LIBRARY 



